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Abstract 

A strong azimuthal flow signature at RHIC suggests rapid system equilibration leading to an 
almost perfect fluid state. The longitudinal extent of the flow behavior depends on how this state 
is formed and can be studied by measuring the pseudorapidity and transverse momentum depen- 
dence of the second Fourier component (v2(pr)) of the azimuthal angular distribution. We report 
on a measurement of identified-particle x>2 as a function of pj (0.5-2.0 GeV/c), centrality (0- 
25%, 25-50%), and pseudorapidity (0 < J] < 3.2) for yfs^ = 200 GeV Au+Au collisions. The 
BRAHMS spectrometers are used for particle identification (n, K, p) and momentum determina- 
tion and the BRAHMS global detectors are used to determine the corresponding reaction-plane 
angles. The results are discussed in terms of the pseudorapidity dependence of constituent quark 
scaling and in terms of models that develop the complete (azimuthal and radial) hydrodynamic 
aspects of the forward dynamics at RHIC. 



In a non-central collision of two relativistic heavy-ions, strong pressure gradients are set up 
in the almond-shaped interaction region. Early results from RHIC have shown that this overlap 
region behaves as an almost perfect fluid, with the greater pressure that exists at the waist of the 
almond leading to greater particle production near the reaction plane. This azimuthal asymme- 
try in particle production is characterized by the strength of the second Fourier component (V2) 
of a harmonic expansion of the angular distribution. The outward pressure can also lead to an 
outward radial flow of the streaming particles, giving a velocity boost to these particles. To fully 
understand the hydrodynamic properties of RHIC collisions, it is necessary to determine the in- 
tegral and differential v'2 behavior, as well as establish the particle spectra and relative particle 
yields [ 1 ] . By exploring the interplay of elliptic and radial flow at both mid- and forward pseudo- 
rapidity, the BRAHMS results better constrain models of the initial conditions and longitudinal 
extent of the interaction region. 

The measurement explores the correlation of charged hadrons and identified particles de- 
tected in the two BRAHMS spectrometers with respect to reaction planes deduced using four 
azimuthally symmetric rings of detectors arranged around the beam line as part of the experi- 
ment's multiplicity array |2y. The V2(pr) dependence of particles detected in the spectrometers 
was determined by the standard reaction-plane method [3]. The reaction-plane resolution cor- 
rection was based on a full GEANT simulation of the BRAHMS experimental response. A 
pseudoevent generator was used to obtain a particle throw consistent with previously established 
particle spectra measured using the BRAHMS spectrometers, with the azimuthal asymmetry of 
the particle throw set to reproduce the PHOBOS integral v'2 results [4]. The correction factor was 
taken as the ratio of the V2(pr) values obtained from the reconstructed events with that input to 
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Figure 1: Charge-hadron V2{pr) for central and mid-central events for r\ as 0, 1 and 3. The curves show 3D Hy- 
dro+Cascade calculations |5|] ■ 

the event generator. 

Figure Q] shows the resulting charged-hadron V2(pr) values for 0-25% and 25-50% central 
events. The increasing v>2 behavior up to pj « 1.5 GeV/c is characteristic of hydrodynamic 
flow Near mid-rapidity the results show very little pseudorapidity dependence. In going 
to forward pseudorapidity with r\ « 3, the slope is found to decrease for both centrality selec- 
tions, with a greater decrease for the mid-central events. The curves show the results of 3D Hy- 
dro+Cascade calculations employing Glauber motivated initial conditions 0]. Good agreement 
is found with experiment for all but the mid-central, forward-pseudorapidity results. Turning off 
the hadronic cascade part of the calculation leads to near pseudorapidity independence of the 
vi(Pt) results, in strong disagreement with the observed slopes at rj « 3. Folding the Fig. Q] 
results with the corresponding charged hadron spectra measured at rj » and 3 yield integral 
V2 values of 0.036 ± 0.005 and 0.027 + 0.004, respectively, at the two pseudorapidities, in good 
agreement with the PHOBOS integral V2 results 14j. 

The identified particle vtipr) results are shown in Fig. [2] Again, the 3D Hydro+Cascade 
calculations, as shown by the smooth curves, are in good agreement with experiment except 
for the forward-pseudorapidity, mid-central events. In addition to reproducing the experimental 
slopes, the calculations also do a good job reproducing the observed mass ordering. 

One of the remarkable features of the RHIC elliptic flow behavior is how closely it follows 
that expected for a perfect fluid. This is particularly evident when the elliptic flow V2(pr) values 
are scaled by the eccentricity e of the overlap region and the number of consituent quarks of the 
detected particle n q and then plotted against the mean transverse energy per constituent quark, 
(£Y)/« 9 , as shown in Fig.[3]for the BRAHMS results at rj «0, 1, and 3. With this scaling, the 
elliptic flow observed for a/^nn = 200 GeV Au+Au collisions for a large number of different 
outgoing particle types are found to follow a common trend [7], as shown by the curve in the 
figure. Such behavior is consistent with the creation of a near-perfect fluid, with the constituent 
quark scaling suggesting this fluid involves quark degrees-of -freedom, as expected in coales- 
cence models. Our central results are consistent with the established systematics at all three 
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Figure 3: V2 scaled by the number of valence quarks n q and the participant eccentricity 6 as a function of the transverse 
kinetic energy KEx scaled by n„. The curve, common to all panels, is based on the systematics observed for a large 
number particle types at RHIC 0|. 



pseudorapidities. For mid-central collisions, the data continue to track well with the systematics 
for r] and 1, but show significantly reduced elliptic flow at forward pseudorapidity. As also 
found in comparison of 3D Hydro+Cascade results to our data, the mid-central events at forward 
pseudorapidity suggest a process other than ideal hydrodynamics is playing a role. 

The elliptic flow behavior is believed to be established at an early stage of the reaction, with 
the integral elliptic flow largely fixed at the point of chemical freeze out |H]. The differential 
elliptic flow also depends on the subsequent hadronization stage, where radial flow can signif- 
icantly affect the final particle spectra. Radial flow results in a velocity boost of the outwardly 
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Figure 4: Transverse kinetic energy (Et) as a function of the mass for pions and protons at rapidity y»0 and 3. The 
curves are linear fits to guide the eye. 

streaming particles. Figure [4] shows the observed {Et) values for pions and protons at rapidities 
y=0 and y«3. It is observed that the transverse energy decreases in going to forward rapidity, 
suggesting a possible reduction in the radial flow. This change in the radial flow can strongly 
influence the differential elliptic flow behavior, but is expected to have less of an effect on the 
integral flow. 

In conclusion, BRAHMS has measured identified particle V2(pr) at rj =s0, 1, and 3 for the 
Au+Au system at y/s^ = 200 GeV. The differential elliptic flow decreases at forward pseu- 
dorapidity, with the decrease for central events consistent with the expectations of 3D Hy- 
dro+Cascade calculations. For mid-central collisions at forward pseudorapidity the elliptic flow 
is found to be significantly less than expected by hydrodynamic calculations. This reduction in 
the elliptic flow is also evident when the current results are compared to previous mid-rapidity 
results using constituent quark scaling. A decrease is observed in the mean transverse energy of 
particles going to forward pseudorapidity, suggesting a reduction in the radial flow component. 
This change in radial flow has a significant influence on the differential V2(pr) values. In general, 
the forward pseudorapidity elliptic and radial flow results place significant constraints on rapidity 
dependent model calculations of the dynamics of RHIC collisions. 
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